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ABSTRACT 

We present a summary of recent assessments of the mass distribution in disk galaxies. Of issue 
in order to characterise galaxy formation models is to determine the relative fraction of baryons and 
dark matter at all radii in galaxies. For disk galaxies, various measurements of the mass distribution 
in galaxies based on vertical kinematics, strong lensing, residuals of scaling relations, fluid dynamical 
modeling, bar strength and pattern speed, warps, and others, have called for either a maximal or 
sub-maximal contribution of the baryons in the inner parts of the disk. We propose a global picture 
whereby all galaxies are typically baryon-dominated (maximal) at the center and dark-matter dom¬ 
inated (sub-maximal) in their outskirts. Using as a fiducial radius the peak of the rotation curve of 
a pure baryonic exponential disk at R2.2 = 2.2 Rd, where Rd is the disk scale length, the transition 
from maximal to sub-maximal baryons occurs within or near R2.2 for low-mass disk galaxies (with 
Wot ~ 200 kms -1 ) and beyond R2.2 for more massive systems. The mass fraction’s inverse depen¬ 
dence on circular velocity in disk galaxies is largely insensitive to the presence of a bar at I? 2 . 2 - The 
mean mass fractions of late- and early-type galaxies are shown to be different at the same fiducial 
radius and circular velocity, pointing to different galaxy formation mechanisms. Feedback, dynamical 
friction, size, and IMF variations are key ingredients for understanding these differences. 

Subject headings: galaxies: spiral - galaxies: elliptical - galaxies: formation - dark matter 


1. INTRODUCTION 

Ever since the first spec ulations about the presen ce of 
“dark matter” in galaxies (lKaptevnlll922HOortlll932lf . as¬ 
tronomers have endeavored to characterize the fraction 
of visible to total mass at all radii in galaxies. The rela¬ 
tive distribution of matter in galaxies ought to be one of 
the most definitive predictions of galaxy formation mod¬ 
els yet its validation is challenged by numerous obser¬ 
vational, theoretical, and operational challenges. While 
galaxies are believed to be dominated by an invisible 
matter component in their outskirts, a debate has been 
blazing for the last two decades regarding the relative 
fraction of baryons and dark matter in the inner parts of 
galaxies: whether galaxies are centrally dominated by 
baryons (“maximal disk”) is of issue. Some of those 
debates have been misconstrued on account of opera¬ 
tional confusion, such as dark matter fractions being 
measured and compared at different radii. An ultimate 
goal for galaxy structure studies is to achieve accurate 
data-model comparisons for the relative fractions of bary¬ 
onic to total matter at any radius. 

In this Letter, we wish to provide a global apprecia¬ 
tion of the mass distribution of baryons and dark mat¬ 
ter in disk galaxies by combining various reports of the 
baryon velocity fraction, R(R) = War/Wot in disk galax¬ 
ies, where War is the inferred circular velocity ascribed 
to the baryonic material and Wot is the total measured 
circular velocity at a given radius. W>ar is typically in¬ 
ferred by multiplying a galaxy light profile, L(R), by a 
suitable mass-to-light ratio, (M*/L) ( R ). Any difference 
between War and Wot is viewed as a signature for the 
presence of dark matter. The dark matter fraction is de- 
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fined as /dm(-R) = / M tot (R) ~ 1 --A 2 , where 

Mdm(R) and M tot (R) are the dark matter and total mass 
enclosed within radius R , respectively. The approxima¬ 
tion becomes an equality for spherical systems. 

How galaxies ultimately arrange their baryonic and 
non-baryonic matter is the result of numerous complex 
mechanisms involving their mass accretion history, the 
depth of the potential well, the initial mass function 
(which ultimately affects feedback and quenching pro¬ 
cesses), dynamical friction, dynamical instabilities, and 
more. Analytical models of galaxy form a tion strive to 
predict fnJ(R) (iDalcanton et al. 1 Il997t iDutton et al.l 
l2007i : Mo et all 120101 1 with mitigated success in light of 
the above challenges; numerical simulations suffer addi¬ 
tional limitations too. 

Can observations provide unique determinations of 
/dm(-R) as ideal constraints to galaxy formation mod¬ 
els? Not yet, but, after decades of muddled debates and 
incomplete data, a clear picture seems to be emerging. 

2. MASS FRACTION ASSESSMENT 

Mass models of disk galaxies, where the observed 
rotation curve is decomposed into its principal gas, 
stars, and dark matter components, have historically 
been embraced as the ideal method to sep arate baryons 
and dark matter as a function of radius dBo sma 1978; 

iCarignan fe Freemanlll985tlvan Albada et al.lll985li . The 

pros and cons of this approa ch are addressed in th e re¬ 
view on “Galaxy Masses” by ICourteau et al. £201 il l. As 
stressed in that review (see also lDutton et al.ll2005H . fun¬ 
damental degeneracies between the disk and dark halo 
models prevent a unique baryon/dark matter decompo¬ 
sition based on rotation curve data alone. Other infor¬ 
mation or methods must be considered to break these 
degeneracies. 

Among others, the quantity War is obtained indirectly 
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via a stellar or baryonic mass-to-light ratio (M*/L) which 
is itself inferred from stellar population studies or dynam¬ 
ical stability analysis. Current stellar population rnass- 
to-light ratios carry an uncertainty of a factor 2 at best 
(IConrovI 120131: ICourteau et al.l 12014ft . The latter uncer¬ 
tainty is unfortunately large enough to encompass a full 
range of baryon fractions at a given radius. 

For reference, a maximal disk is often defined to obey 

F = Vdisk {R max )/V tot (i? max ) > 0.85, (1) 

where Vdisk is the inferred velocity of the disk (stars 
and gas), Vtot is the total observed velocity, and i? max 
is the radius at which lAisk reache s its peak value 
(|van Albada et al.1l 198.51 ISackettiri997l) . 

In terms of dark matter fraction this criterion corre¬ 
sponds to /dm < 0.28. 

2.1. Tully-Fisher residuals 

To circumven t the loopholes of mass modeling, 
ICourteau fe Rbd (I1999T) observed that the residuals of 
the velocity-luminosity relation, AlogV(L), contrasted 
against those of the size-luminosity relation, A log R(L), 
can provide a tighter constraints on T measured at 
2.2 i? d . 

It can be easily shown that A log V ( L)/Al o gR(L) = 
— 1/2 for baryon dominated galaxies (ICourteau k, Rbd 
Il999h and A log V(L)/A log R{L) > 0 for dark 

matter dominated galaxies. Numerous scaling re¬ 
lation analyses of spiral gala xies have found that 
A log V(L) /A log R(L ) cs 0 (ICourteau fc R.ixl IT999I: 
Courteau et al.l 120071 : iDutton et al.l 120071 iReves et al.1 
20111) . The null slope indicates that both baryons and 
dark matter contribute to the galaxy dynamics at 2.2 Rd- 
Various galaxy structure models whereby the dark mat¬ 
ter is compressed adiabatically by the cooling baryons 
as the galaxy stabilizes dynamically yield a mapping be¬ 
tween A log V (L) / A\og R(L) and the baryon fraction T 
(see references above). These models yield, on aver¬ 
age for bright spiral galaxies, T(2.2 Rd) = 0.55 — 0.71. 
The lower value applies if only adiabatic contraction 
is at play; a higher value is found if other mecha¬ 
nisms (e.g. feedback) oppose the adiabat ic contraction 
(|Courteau fe Rbd[l999HDutton et al.ll2007D . By virtue of 
Eq. m above, typical galaxy disks are thus, on average , 
sub-maximal at 2.2 Rd- Since mass scale as V 2 , t he d ark 
matte r fraction, /DM(2.2i?d) — 50 — 70%. IReves et al.l 
(|2011f) suggested that T is controlled by the surface den¬ 
sity of the baryonic material; i.^ if the latter is high, so 
is J 7 (see also Fig. 23 in lBovv fc Rbdl2013 /. 

2.2. Vertical Disk Kinematics 

A similar result is ob tai ned via the “velocity disper¬ 
sion” method of (Ivan der Kruitl 11988 / . as app li ed b y 
many for a full s uite of disk galaxie s (IBottemal 119931) . 
edge-on systems dKregel et al.l 20051) . or face-on spira l 
galaxies dBershadv et al.l 120111 : Martinsson et al.ll2013D . 
For a self-gravitating, radially exponential disk with ver¬ 
tical profile o f the f orm p{R,z) = p(R, 0)sech 2 (z/z n ), 
Ivan der Kruitl (|1988l) and IBottemal (1199311 showed that 
the peak circular velocity of the stellar disk, Vusk, mea¬ 
sured at R = 2.2 can be related to the vertical ve¬ 
locity dispersion, V z , and the intrinsic thickness (or scale 


height) of the disk, z o, via: 

Vdi S k(Rpeak) = C (V z 2 )^ 0 (2) 

V zo 

where c ~ 0.88(1 — 0.28zo/Rd) dBershadv et al.l 1201 ID . 

Edisk(7?peak)j the contribution of the baryonic disk to 
the total velocity where the rotation curve reaches its 
peak (at i? pea k), can thus be measured directly if the 
disk scale length, l?d, the scale height, zq, and the verti¬ 
cal component of the velocity dispersion, V z , are known. 
For face-on systems, V z and Rd can be measured but zq 
must be inferred; and vice versa for edge-ons. Van der 
Kruit’s method is thus statistical in nature and relies on 
a n um ber of (p oten tially noisy) scaling transformations. 
iBovv fe Rbd (12013H also point out that the velocity dis¬ 
persions measured from integrated light, as in most ap¬ 
plications of Eq. (J2J for external galaxies so far, do not 
trace the older, dynamically-relaxed stellar populations 
with high fidelity. The younger, brighter, stellar popula¬ 
tions traced by the light-weighted spectra tend to have 
a lower velocity dispersion therefore biasing Vdisk(Rpe a k) 
low and thus favoring sub-maximal disks. Sample selec¬ 
tion also plays a substantial role (see Fig. 1 below). 

The DiskMass project has made use of van der 
Kruit’s formalism to extract Vdisk (Rpea,k) for a sample 
of 46 nearly face-on (inclinations ~ 30 degrees) galax¬ 
ies wit h rotation velocities b etween 100 kms -1 and 250 
kms~ 1 (|Bershadv et al.ll2011 ,). This survey uses integral- 
held spectroscopy to measure stellar and gas kinemat¬ 
ics using the custom-built SparsePak and PPAK in¬ 
struments. For 30 Disk Mass galaxies covering a range 
of structural properties, iMartinsson et al.l (|2013l ) report 
that the fraction T ranges from 0.34 to 0.73 and increases 
with luminosity, surface brightness, rotation speed, and 
redder color. The average val ue of their sample is T = 
0.57 ± 0.07, in agreement with ICourteau h Rixl (I1999D . 

2.3. Other methods 

Similar and broader results have been achieved through 
very different techniques for disk domina ted systems : 
such as fluid dynamical modeling dWeiner et al.l 120011 : 
iKranz et alJl2003l : IAthanassoulall2014l see the latter ref¬ 
erence for a more critical outlook) and strong gravita - 
tional lensing (lBarnabe et al.ll20l2 IDutton et alJl2013aj) . 
For instance, because gravitational lensing is sensitive to 
projected mass while dynamics trace the mass within el¬ 
lipsoids, their combination provides additional geomet¬ 
rical dimensions to disentangle the baryons and dark 
matter in galaxies; this approach is especially efficien t 
for highly-inclined disk gal axies (IDutton et al.l 1201 lbl) . 
These methods, reviewed in ICourteau et al.1 (120111 ). are 
consistent on average with ,F(2.2 Rd) = 0.6 —0.7, but also 
show a broad range of T as a function of maximum circu¬ 
lar velocity or disk size. In their hydrodynamical model¬ 
ing of five grand design non-barred galaxies. IKranz et al.1 
( 2003 / presciently suggested that galaxy disks appear to 
be maximal if Vn a x > 200 kms -1 , sub-maximal other¬ 
wise. The later summarizes rather well the more com¬ 
plete picture about the mass distribution in galaxies that 
is now clearly emerging. 

3. A UNIFIED VIEW 
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Fig. 1.— Dark matter fraction at 2.2 as a function of total 
circular velocity for samples discussed in the text. A clear trend is 
observed such that the dark matter fraction decreases with circular 
speed; indeed below 250 kms -1 , dark matter dominates (> 50%) 
the total mass budget of th e galaxy at 2.2 R< -\ . The yellow shading 
shows the stability criteria dEfstathiou et alTl982D for stellar disks 
(/dm — 0.67; short-dashed black line and above) and gaseous disks 
(/dm — 0.51). The dotted black line defines a maximal disk (Eq. 
[1]). The black line represents a model for g alaxy scaling relation s 
assuming a Chabrier initial mass function dDutton et al.|[2011aO ; 
the long-dashed black lines are the lcr dispersion in galaxy sizes. 


Fig. |T| is a compilation of dark matter fractions, /dm = 
1 — T 2 (W>M/Wot) 2 , at ~ 2.2 i?d for various galaxy 
samples discussed above (others are discussed below), as 
a function of total circular velocity. 

The DiskMass survey (blue points) predominantly 
samples large-to-average size galaxies thus avoiding the 
regime of disk maximality (at 2.2 i?d)- The cyan pen¬ 
tagons represent the pa rameter space for the T ully-Fisher 
residual analysis of the iCourteau et al.l (120071) sample of 
spiral galaxies with (open pentagon) and without (filled 
pentagon) adiabatic contraction (IDutton et al.1 120071) . 
The Courteau and DiskMass samples target essentially 
the same sub-maximal galaxy types. Unlike those sam¬ 
ples, the SWELLS gravitational lensing s urvey is intrin- 
sically biased towards high mass systems (Barnab e et al.1 
2012t IDutton et aTI l2013al) . The Milky Way study of 
Bow & Rix (120131 ) uses a method similar to that of 
van der Kruitl (119881 ) but with the added benefit that 


stellar populations can be disentangled, thus providing 
less biased measurements of the true underlying structure 
(scale length, rotation and dispersion) of the dynamical 
tracer. Their analysis results in a low mass-weighted disk 
scale length, = 2.15±0.14 kpc, and a correspondingly 
low estimate of /dm = 0.31 ± 0.07 at 2.2 i?J3- 
The short-dash boundary delineates the region 


3 Various investigation s of the M ilky Way m ass density profiles 
ar e reported in Table 1 of Courteau et al. (2014). The ea r ly stu dies 
of IDehnen fc Binnevl 119981 ') and [Englmai er’&r Gerhardl 11999H al¬ 
ready pointed to a maximal MW disk. Measurements of T for the 
Milky Way hinge significantly on estimates of its disk scale length. 
A low estimate of R r { yields a high value of T (or /dm)- 


where baryon-heavy disk galaxies become disk un¬ 
stable. The criterion for stability against a bar 
is Vc.wci Rd /GM^ ) 1 / 2 £ 1.1 dEfstathiou et al.l 119821 : 
iFovle et all I2008D . The yellow shading shows the sta¬ 
bility criteria for stellar (/dm > 0.67) and gaseous 
(/dm > 0.51) disks. Barred and unbarred galaxies have 
also been shown to have the same velocity-luminosity 
and size-luminosity sc a ling relations (mea sured at R 2 . 2 ) 
(|Courteau et al.l 120031 : iSheth et al.1120121) . Thus, in a 
global dynamical sense, barred and unbarred galaxies 
behave similarly and have, on average, comparable frac¬ 
tions of luminous and dark matter at a given radius. 

The black lines in Fig. [T] show a toy model w h ere we 
have used the observed sca ling relations dCourteau et al.1 
120071 : IDutton et al.ll20lTal) between disk scale length, ro¬ 
tation velocity and stellar mass (assuming a Chabrier 
IMF) to estimate the dark matter fraction, assuming the 
baryonic mass is dominated by an exponential stellar disk 
(i.e., gas and stellar bulges are ignored). The dashed lines 
show the effect of changing the disk sizes by ±0.16 dex 
(the scatter in the size-velocity relation). Changing the 
stellar masses (e.g., by changing the IMF) by the same 
amounts leads to similar changes in dark matter frac¬ 
tions. This simple model captures the main features of 
the observations: 

1. Dark matter fraction decreases with galaxy circular 
velocity; 

2. Milky Way mass galaxies (Wire ~ 220 kms -1 ) 
have, on average at 2.2 /dm ~ 0.5; and 

3. There is significant scatter in dark matter fractions 
at a given circular velocity driven, at least partly, 
by size variations. 

3.1. Early-Type Galaxies 

Measuring dark matter fractions in early-type galax¬ 
ies is more challenging. Firstly, because measuring total 
mass profiles is subject to the mass-anisotropy degener¬ 
acy, which requires high-quality data to break. Secondly, 
because there are no direct dynamical means of separat¬ 
ing the total mass into baryons and dark matter. Ob¬ 
servations of dark matter fractions in early-type galaxies 
are thus driven by systematics in the mass modeling as¬ 
sumptions, such as dark matter profile, and stellar M/L 
gradients. 

Fig. [H also shows the distri bution /pm(-Rb) of 26 0 
ATLAS3D early-type galaxies (|Cappellari et al.112013D . 
where we have used Wire = 1-38(78 to tran sform veloci ty 
dis per sions, cr 8 , to circular velocity Wire (IDutton et all 
l2011bl) . An important caveat to these results is the as¬ 
sumption of an NFW dark matter halo (Navarro, Frenk, 
& White 1997) in the mass modelling. Dark matter frac¬ 
tions as high as /dm (-Re) ~ 60% are obtained if dark mat¬ 
ter haloes contract adiabatical ly in response to galaxy 
formation (IDutton et al.l l2013bl) . Note also that AT- 
LAS3D velocity maps do not typically extend beyond 1- 
1.3 R e , the regime reported for LTG in measuring V 2 . 2 - 
Therefore, the data for ETGs presented in Fig. |T| are 
likely lower limits; still, the latter probe ETG dynamics 
in a regime of maximal ‘baryons”. 

Fig.[2]shows dark matter fractions vs radius for ACDM 
based models of early and late-type galaxies with circu- 
lar velocity Wirc(-Re) — 230kms -1 from IDutton et all 
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R/Re 


Fig. 2.— Dark matter fraction as a function of half-light radius, 
R e , for early-type (red lines) and late-type (blue lines) galaxies with 
circular velocity V^i rc (i? e ) — 230 kms -1 . The black horizontal lines 
are as in Figure 1. Solid lines show models with NFW ha loes that 
follow the concentration-mass relation from ACDM (1 Maccio et al.l 
120081) . Dashed lines show models with adiabatically contracted 
NFW haloes and correspondingly lighter IMFs. The blue dots 
correspond to 2.2 R d for late-type galaxies. 

( 2011a ). This shows that for a given halo response, 
early-type galaxies have lower dark matter fractions than 


late-type galaxies. These differences are driven by both 
size (early-types are smaller) and IMF (early-types have 
heavier IMFs) variations. 

4. FUTURE OUTLOOK 

The trends observed in Fig. |T] and the characteriza¬ 
tion of baryon-DM transitions in galaxies over a range of 
Hubble types are significant constraints that galaxy for¬ 
mation models have yet to uniquely predict (e.g. Fig. 0. 
This is partly due to additional complications such as the 
degeneracy between the baryon or dark matter fractions, 
the stellar initial mass function, adiabatic contraction, 
and various counter-acting non-gravitational agents suc h 
as feedback and dynamical friction dD utton et al.ll20liat 
iTruiillo-Gomez et al.l 120111 : iQguri et al.112014T ) . The ad¬ 
vent of new extensive dynamical models of galaxies based 
on wide-held integral held surveys, such as CALIFA, 
MaNGA, SAMI (out to ~ 2 R e ) and SLUGGS (reach¬ 
ing ~ 6 — 10f? e ) heralds a promising future for mapping 
the mass distributions at all radii in galaxies. Model de¬ 
pendencies and biases could at least be partially offset 
with reduced sampling errors. 
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